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Abstract 

We examine the spinless one-dimensional Falicov-Kimball model (FKM) below half-filling, ad- 
dressing both the binary alloy and valence transition interpretations of the model. Using a non- 
perturbative technique, we derive an effective Hamiltonian for the occupation of the localized or- 
bitals, providing a comprehensive description of charge order in the FKM. In particular, we uncover 
the contradictory ordering roles of the forward-scattering and backscattering itinerant electrons: 
the latter are responsible for the crystalline phases, while the former produces the phase separation. 
We find an Ising model describes the transition between the phase separated state and the crys- 
talline phases; for weak-coupling we present the critical line equation, finding excellent agreement 
with numerical results. We consider several extensions of the FKM that preserve the classical na- 
ture of the localized states. We also investigate a parallel between the FKM and the Kondo lattice 
model, suggesting a close relationship based upon the similar orthogonality catastrophe physics of 
the associated single-impurity models. 

PACS numbers: 71.10.Fd, 71.30.+h 



1 



I. INTRODUCTION 

The Falicov-Kimball Model (FKM) describes the interaction between conduction elec- 
trons and localized atomic orbitals. The Hamiltonian of the one-dimensional (ID) FKM for 
spinless Fermions is written 



where t > is the conduction (c) electron hopping, e/ is the energy of the localized /-electron 
level, and G is the on-site interorbital Coulomb repulsion. The concentration of electrons 



consider only the case n < 1. We work throughout at zero temperature T = 0. 

The FKM was originally developed as a minimal model of valence transitions: continuous 
or discontinuous changes in the occupation of the / orbitals (the atomic "valence") were 
observed when varying the coupling G or the /-level energy e/.— Since only the distribution 
of electrons across the two orbitals is of interest, the model has traditionally been studied for 
spinless fermions. These early works, however, neglected an important feature of T^fkm- the 
occupation of each /-orbital is a good quantum number and so may be replaced in Eq. by 
its expectation value rij —>■ {rij) = 0, 1. It was quickly realized that in many physical systems 
displaying a valence instability (e.g. SmBg and Ce), this is an inappropriate idealization. 
Instead of a mixture of atoms with different integer valence, in these materials each atomic 
orbital exists in a superposition of its different occupancy states.— Although the FKM was 
modified to include this quantum behaviour by the addition of a c-f hybridization term,^^ it 
has now been superseded as a model of valence transitions by the periodic Anderson model.— 

The FKM was reinvented by Kennedy and Lieb in 1986 as a simple model of a binary 
alloy.— Assuming fixed c- and /-electron populations, the sites with occupied and unoccupied 
/ orbitals may be regarded as different atomic species A and B respectively. The Coulomb 
repulsion G is interpreted as the difference between the single-particle energies of the two 
atoms. For this so-called crystallization problem (CP), the ground state is defined as the 
configuration of the two atomic species (/ electrons) that minimizes the energy of the c 
electrons. The ordering of the different atomic constituents in a binary alloy is an important 
theoretical and experimental problem: in realistic systems a large range of ordered structures 
are observed, although the electronic mechanisms responsible for these phases have remained 




(1) 




where N is the number of sites. In this work we 
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largely obscure.i^ By studying a simple model such as the FKM, some insight into the origin 
of the charge order might be obtained. 

Kennedy and Lieb analyzed Eq. (^Q) for a bipartite lattice at half-filling and equal con- 
centrations of c and / electrons. In the limit of T = and strong-coupling, they proved 
that the / electrons occupied one sublattice only, the so-called checkerboard state. This 
crystalline state is, however, unique to half-filling: for all other fillings, the G — oo ground 
state is the so-called segregated (SEG) phase.— The SEG phase is characterized by the 
/ electrons forming a single cluster, arranged in such a manner as to present the smallest 
perimeter with the rest of the lattice, which is occupied by the c electrons. These strong- 
coupling results hold for all dimensions d. At weak- and intermediate-coupling, the situation 
is considerably more complicated: for d = 1, both analytio^i^ and numeri o^^i^^i^^ studies 
have revealed a myriad of different crystalline orderings of the / electrons. The SEG phase 
is also realized, but not as ubiquitously as at strong-coupling. Intriguingly, for certain c- and 
/-electron fillings, the system is unstable towards a special phase-separated state, where the 
ground state configuration of the localized electrons is a mixture of a crystalline phase and 
the state with completely empty or full localized orbit als."^^*^^*^^ Work in higher dimensions 
has revealed similar behaviour;^® the understanding of the — > oo limit phase diagram is 
particularly advanced.— 

Contemporary with Kennedy and Lieb's work, Brandt and Schmidt introduced the FKM 
as an exactly-solvable model of a "classical" valence transition.— The distribution of the 
electron weight between the two orbitals is not fixed, but instead determined by the inter- 
ations. Quantum effects such as superposition of orbital states are ignored: as in the CP, 
the valence transition problem (VTP) is also concerned with the configuration adopted by 
the available / electrons. Despite the similarity between the two interpretations, apart from 
the d —* oo limilfi^ and the d = 1 half-filling case ,^^i^^ very little is known about the ground 
states of the VTP. Since the ordered configurations found for the CP occur over a wide range 
of different fillings and coupling strengths, we can nevertheless expect that the VTP has a 
similarly rich phase diagram. 

Although an impressive catalog of charge-ordered phases has been assembled for the 
ID FKM, only the weak-coupling crystalline phases are easily explicable as due to the c- 
electron backscattering off the localized orbitals. The mechanism responsible for the weak- 
coupling segregated and phase separated states remains unknown; the competition between 
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crystallization and the segregation is also poorly comprehended. In this paper, we expand 
upon our previous work,— outlining a comprehensive theory of the charge order in the FKM, 
with particular emphasis on the phase separation. 

To describe the c electrons, we use the well-known non-perturbative bosonization tech- 
nique, specially adapted to account for the presence of localized orbitals. We then canonically 
transform the bosonized FKM, rewriting the Hamiltonian in a new basis that reveals the 
origin of the phase separation to be the c-electron delocalization. Such a mechanism has 
previously been proposed to account for the ferromagnetic phase in the ID Kondo lattice 
model (KLM),— pointing to a nontrivial connection between the two models based upon 
orthogonality-catastrophe physics. After simple manipulation of the transformed Hamilto- 
nian we decouple entirely the c and the / electrons, obtaining an Ising-like effective Hamil- 
tonian describing only the occupation of the / orbitals. The competition between the seg- 
regation and crystallization is clearly evident in this effective model: at weak-coupling we 
find the backscattering crystallization dominates the physics; with increasing G, however, 
the electron delocalization drives the system into the SEG phase. We verify that both 
crystallization and segregation are present also in the VTP. 

Our paper is arranged as follows: in Sec. |H] we give a brief outline of our bosoniza- 
tion procedure, and present T^fkm in the bosonic form. We proceed to a description of the 
canonical transform in Sec. IIIIl including a discussion of the resulting terms. We argue 
in Sec. IIVI for the derivation of the effective Hamiltonian for the localized / orbitals from 
the canonically-transformed Hamiltonian; this is subsequently used in Sec. El to interpret 
the numerically-determined phase diagrams for the CP (Sec. IV Ajl and the VTP (Sec. IV B| . 
We also present a brief analysis of several extensions of the FKM in Sec. IVH specifically 
intraorbital nearest-neighbour interactions and the introduction of spin, focusing upon pos- 
sible alteration of the CP phase diagram. We conclude in Sec. IVIII with a summary of our 
results and the outlook for further work. 

II. BOSONIZATION 

The technique of bosonization has for many years been used to study the critical prop- 
erties of one-dimensional many-electron systems.-^ It relies upon the remarkable fact that 
an effective low-energy description of such systems may be constructed in terms of bosonic 
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fields: this representation is usually much easier to manipulate than the equivalent fermionic 
form. The bosonization of a tight-binding Hamiltonian is often performed in the continuum 
limit where the lattice spacing a — >^ 0;— this approach is, however, inappropriate for systems 
involving localized electron states. For the itinerant electrons in the FKM, however, the 
usual bosonization approach can be generalized to account for the presence of the local- 
ized / electrons. As explained in Ref. this is accomplished by imposing a finite cut-off 
a > a on the wavelength of the bosonic density fluctuations. Below we summarise our 
methodology. 

The Bose representation is most conveniently written in terms of the dual Bose flelds. 
For a system of length L ^ a we have 

= -^J2Il ^PMK{k)e''''' (2) 

At the core of the bosonization technique are the chiral density operators 

Pu{k) = q Yl ^k'-kCk' (4) 

0<uk'<-K/a 

which describe coherent particle-hole excitations about the right and left Fermi points: as 
subscript (otherwise) we have v = R{+), L{—) respectively. The Pu{k) are the basic bosonic 
objects, obeying the standard commutation relations 

ukL 

[Pu{k),p^'{k')]_ = 6yy5k-k'^r~ (^) 

ZTT 

for wave vectors < ^. The physical significance of the Bose fields is as potentials: dx4>{xj) 
and dxO{xj) are respectively proportional to the departure from the noninteracting values 
of the average electron density and current at Xj. 

The bosonic wavelength cut-off is enforced in Eq. Q and Eq. (jS)) by the function A^{k) 
which has the approximate form 

Kik) ^ ' ' " (6) 
I otherwise 

We expect that A^{k) is a smoothly varying function of \k\, reflecting the gradual change 
in the nature of the density fluctuations. We require, however, that A^(/c) be not too 'soft', 



i.e. ^^{k) ~ Ao(^) for m = 2,3,4. The cut-off essentially 'smears' the Bose fields over 
the length a below which the density operators do not display bosonic characteristics. The 
commutators of the Bose fields reflect this smearing, with important consequences for our 
analysis: 

[(t){Xj),e{x.y)]_ = —Sgn^{Xj>-Xj) (7) 

[d^(f){xj),9{xj>)]_ = -m5a{xji - Xj) (8) 

sgnQ,(x) and 5a{x) are the a-smeared sign and Dirac delta functions respectively. The precise 
form of these functions depends upon Aq,(/c) [see Sec. 1111 Aj . 

As is customary, we linearize the c-electron dispersion about the two Fermi points. This 
allows a decomposition of the j-site annihilation operator in terms of states in the vicinity 
oi kp (the right-moving fields) and —kp (the left-moving fields): 

Remarkably, the density operators Pv{k) generate the entire state space of the linearized 
Fermion Hamiltonian. A Bose representation for the c^j may then be derived by requiring 
that it correctly reproduces the Fermion anticommutators and noninteracting expectation 
values. This leads to the fundamental bosonization identity 



Aa - 

Cuj = Y exp {-iiy [4>{xj) - Jy6{xj)]) (9) 

We note that this identity is only rigorously true in the long- wavelength limit; Eq. Q may 
not correctly reproduce the short-range (< 0{a)) properties of the c^j- The dimensionless 
parameter A is a normalization constant dependent upon the cut-off function. The Klein 
factors obey the simple algebra: 



Fu^ Fu' 



2F,FA 



F F' 



25,y (10) 



The Klein factors act as "ladder operators": since the Bose fields only operate within sub- 
spaces of constant particle number we require operators to move between these different 
subspaces if we are to regard equation Eq. as an operator identity. That is, -f^ may be 
thought of as lowering the total number of //-moving electrons by one. 

Using standard field-theory methods, the bosonization identity Eq. Q may be used to 
derive the Bose representation for any string of Fermion operators. Of particular note is the 



representation for the j'-site occupancy operator, n^: 

TT a ^ — ^ 

V 

The first term on the RHS, ?7,q, is the noninteracting c-electron concentration; the second term 
gives the departure from this value in the interacting system and is due entirely to forward 
scattering {v ^ v processes); the third term is the first order backscattering (z/ —v 
processes) correction. Higher order backscattering corrections are neglected. 

A. The Hamiltonian in Boson Form 

We bosonize the FKM Hamiltonian using the above methodology. Since only the itinerant 
c electrons can be bosonized, we require that there be a finite population in the noninter- 
acting c-electron band. For the CP this simply requires us to assume finite concentrations 
of the two species, "uF and for the c and / electrons respectively, which do not change 
with the addition of the interaction term. For the VTP, we impose the condition that the 
/-level coincides with the Fermi energy in the noninteracting system: as we consider only 
the case n < 1, we limit ourselves to —It < e/ < — 2t cos(7m/a). For e/ outside this range, 
our bosonization approach does not work. We discuss this in more detail in Sec. IV Bl 

Before bosonizing the Coulomb interaction, we re- write the /-electron occupation in terms 
of pseudospin-i operators, — \ = t/ • In the pseudospin representation, spin-| at site j 
indicates an occupied /-orbital and vice versa. For the CP, the condition of constant /- 
electron concentration then translates into a fixed pseudospin magnetization = — 
The use of the pseudospins will considerably simplify the subsequent manipulations. We 
re-write the Coulomb interaction 

j j j 

We have used the requirement of constant total electron concentration to obtain the second 
term. Substituting Eq. (Illj) into Eq. (fT^ . we obtain the bosonized form of the FKM 
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Hamiltonian 



•FKM — 



vpa 



Ga 



(13) 



TT 




For c-electron concentration rf, the Fermi velocity is defined vp = —2tasm{kpa) where 
kp = im'^/a. Note that the parameter only enters into Eq. ()13|) indirectly through vp 
and kp. Since the Klein factor products in the backscattering corrections [the last term 
in Eq. (fT!?|) ] commute with the Hamiltonian, we replace them by their expectation value, 



III. THE CANONICAL TRANSFORM 

The work on the CP has established that the c electrons mediate interactions between 
the / electrons via the interorbital Coulomb repulsion. Very little, however, is known about 
the character of these interactions: here we seek to reveal the electronic origins of the charge 
order by rotating the Hilbert space basis to decouple the c and / electrons. We apply a 
lattice generalization of the canonical transform used by Schotte and Schotte in the X-ray 
edge problem (XEP);^ 



The canonical transform bears a close resemblance to the transform used by Honner and 
Gulacsi in their analysis of the KLM.— This resemblance is not coincidental, but instead 
points to a fundamental similarity between the FKM and KLM which we explore below. 

A major advantage of the Bose representation is that the transformation of the Bose 
operators under Eq. (fT^ may be calculated exactly using the Baker-Hausdorff formula. 
This allows us to carry the canonical transform of the Hamiltonian through to all orders. 
The effect of the transform may be summarized as follows: 




(14) 




(15) 



7) m ' ^ ■' 




(16) 
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U^p{xj)U = p{xj) - I ) 6a{xj - Xj') (17) 



All other operators in Eq. ()13p are unchanged by the transform. In particular, we note that 
the transform preserves the /-configuration, i.e. UWjU = rj. 

The transformation of the derivative of the 0-field [Eq. is of special note, as it makes 
explicit the dependence of the c-electron density p{xj) = — ^dx(j){xj) at site j upon the 
local /-electron occupation: 

J' 

As expected, the effect of the Coulomb interaction is to enhance (deplete) the c-electron 
density where the / orbitals are empty (occupied). As we explain below, this is the origin 
of the observed segregation in the CP. 

Substituting the transformed Bose fields [Eq. (fT3j) and Eq. (fTBj) ] into Eq. (fT^. we obtain 

U^H^^mU = ^ E {(9Mx,)f + {dj{x,)f} + G (nS - i) 5^r/ 

j 3 



^T^'6a{Xj 



Xj>)Tj, 



3,3 

^2GAa ^ ^2 [0(3; .) _ /c(j) - kpx,]) (18) 



a 

3 

where we have introduced the simplifying notation 

^(j) = ^ ^/'Sgna(a;j' - Xj) (19) 

for the string operator in Eq. (|15p. Since the canonical transformation of T^pkm has been 
carried out exactly, it follows that Eq. ()18p is identical to Eq. ()13p. The result of our trans- 
formation is to have re-written TYfkm in a new basis that includes the effective interactions 
between the / electrons. The rest of this paper will be concerned with the study of Eq. ()18p: 
we begin by examining the origins of the terms involving the / electrons in the transformed 
Hamiltonian. 



A. The Ising interaction 

The removal of the forward-scattering Coulomb interaction by the canonical transform 
introduces an effective interaction between the / electrons: 

G^a" 



2tcvp 

3,3' 



J2T^6^ix,-x,,)r^, (20) 



Unlike other effective interactions, such as the weak-coupling RKKY theory^ or the large-G 
expansion,— Eq. (^01) is non-perturbative. Furthermore, Eq. ()20|) differs from these other 
effective interactions in being responsible only for the segregation and phase separation. Its 
significance warrants some discussion upon its properties and origins. 

The interaction is implicitly dependent upon the properties of the c electrons: the form 
of the potential in Eq. (pn|) is the Fourier transform of the cut-off function 



ux,) = ^j2^^(^y'''' (21) 

k 

To concretely illustrate the variation of the interaction, we consider two choices of cut-off 



0(1^1 — -) step function 
Kik) = { (22) 
exp ( — a I /c I ) exponential 

Here 0(x) denotes the well-known Heaviside step function. For simplicity, we evaluate the 
summation Eq. !^21\i in the thermodynamic limit L — oo and for a continuum system (valid 
for a ^ a). We thus find 

f°° dk 

5a{x) = I — cos(A;x)A^(A;) 



'0 



sin(7rx/a;)/(7rx) step function 

(23) 

7r~^a/ (a^ + x^) exponential 

These integrals are plotted in Fig. this plot makes it clear that a characterizes the 
range of the interaction Eq. pO|) . For the step-function cut-off the potential will take neg- 
ative values for x > a. Since a is limited below by the lattice constant, however, the 
nearest-neighbour value 5a{a) is always non-negative and exceeds in magnitude all other 
values of the potential. In the pseudospin language the interaction Eq. (pUj) is ferromagnetic 
below the bosonic wavelength cut-off. Furthermore, beyond this length scale the potential 
is insignificant. 

The canonical transform reveals that the forward-scattering mediates attractive interac- 
tions between the / electrons; as such, it can account for the observed segregation^ and phase 
separation . ^^1^^ This is not unexpected, as the forward-scattering c electrons transfer small 
crystal momentum (<^ kp) to the / orbitals, thus only interacting with the long-wavelength 
features of the underlying /-electron configuration. 

10 




FIG. 1: The form of the interaction potential Eq. H2U|) for step function [A^{k) = Q{\k\ — vr/a)] 
and exponential [A^(/c) = exp(— q|A;|)] cut-offs. 

To fully understand the physical origin of Eq. ()20|) we must consider the details of the 
bosonization process. Because of the bosonic wavelength cut-off, our treatment can only 
describe density fluctuations over distances > 0{a). The bosonic fields cannot distinguish 
separations less than this distance, hence the smeared canonical field commutators Eq. ((7j) 
and Eq. (jH)). Our description of this system thus assumes that the c electrons are delocahzed 
over a characteristic length ~ a: the a-smeared 5-functions [Fig. (^] may be very crudely 
conceived as the probability density profile of these delocahzed electrons, i.e oc Sa{x). 

The finite spread of the c-electron wavefunctions carries the interorbital Coulomb repulsion 
over several lattice sites [see Eq. (fTTj) ]. directly leading to the segregating interaction Eq. ((201) • 

In the familiar bosonization description of one-component systems such as the Hubbard 
model, making a arbitrarily small does not alter the critical properties of the model; in 
particular, we still obtain the Luttinger liquid fixed point behaviour. In these models a is 
regarded as a short-distance cut-off which defines the minimum length scale in the system 
(usually the average inter-electron separation ~ kp-^), much as the infrared cut-off in field 
theory. To understand the long-wavelength behaviour we need only keep a as a formally 
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finite parameter. Such arguments cannot, however, be made for the c electrons in the FKM, 
where the limiting length scale of the bosonic description is determined by the interactions 
with the localized / electrons: the parameter a therefore enters our bosonic theory as a finite 
but undetermined length. We estimate a by examining the short-range fermionic scattering 
of the c electrons off the / orbitals. 

As is well-known, the configuration adopted by the / electrons in the FKM acts as a 
single-particle potential for the c electrons. That is, the c electrons move in a site-dependent 
potential that takes only two values +G/2 or —G/2, corresponding to occupied and un- 
occupied underlying / orbitals respectively. Below the average inter-electron separation 
{k'^^ = a/nn'^), the c electrons move independently of one another and their motion is there- 
fore described by a single-particle Schrodinger equation. In the limit n'^ —>■ the average 
inter-particle separation is much larger than the lattice constant: it is here acceptable to 
take the continuum limit of the lattice model, yielding a simple form for the Schrodinger 
equation describing the low-energy {E = 0) wavefunctions ipi^)'- 

dli,{x) = Gm{T^)i;{x) (24) 

here m is the bare electron mass. 

The motion of the c electrons across the lattice is analogous to the familiar problem 
of elementary quantum mechanics of a particle in a finite well.-^ For a c-electron moving 
in a region free from / electrons ((t|) = — |), the energy of the c electron exceeds the 
potential and so we find the solutions for iIj{x) to be plane waves. In contrast, the c-electron 
energy is less than the potential in a region of occupied / orbitals ((r^) = we therefore 
find exponentially decaying solutions iIj{x) ~ exp(— a;/C) characterized by the length scale 
C ~ a/I/G. Since we identify a with the finite spread of the delocalized c electrons across 
the orbitals, we conclude that a ~ C- We therefore expect a = h^/tjG where 6 is a constant 
to be determined. 

I. Relationship to the KLM 

The similarity of the canonical transform Eq. (fT^ to that used by Honner and Gulacsi in 
their treatment of the KLM suggests a connection between the two models. This relationship 
is best revealed by considering the single-impurity limit of these lattice models; for the FKM, 
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the associated impurity problem is the XEP. 

As is well known, the sudden appearance of the core hole in the XEP excites an infinite 
number of electron-hole pairs in the conduction band, leading to singular features in the 
X-ray spectrum (the orthogonality catastrophe). Schotte and Schotte recast the problem 
in terms of Tomonaga bosons: the core hole potential directly couples to the boson modes 
of the scattering electrons, and may be removed by a suitable shifting of the oscillator 
frequencies.— Our own canonical transform Eq. (fT^ repeats this procedure across the ID 
lattice. Although the / electrons in the FKM are static, the appearance of the core hole in 
the XEP is equivalent to suddenly turning on the interactions in the FKM. Since we start 
with non- interacting boson fields in Eq. (jl3p . this is a perfect analogy. 

The spin-i Kondo impurity is another classic example of the orthogonality catastrophe, 
although the singular behaviour here arises due to the shifting of the spin-sector boson 
frequencies. In the usual Abelian bosonization approach the boson modes only directly 
couple to the 2;-component of the impurity spin: ignoring the transverse terms the problem 
is identical in form to the XEP. Although these transverse terms somewhat complicate the 
analysis, for special values of (the Toulouse point), it is possible to map the problem to 
the exactly solvable resonant-level model by shifting the c-electron boson frequencies as in 
the XEP.^ For the lattice case this argument may be generalized to arbitrary J^:^^ Honner 
and Gulacsi's canonical transform therefore shifts the KLM's spin-sector boson frequencies 
in precisely the same way as the transform Eq. (|THl shifts the charge-sector boson frequencies 
in the FKM. 

The similarity between the charge-sector physics of the FKM and the spin-sector physics 
of the KLM suggests a parallel between the segregating interaction Eq. and the Kondo 
double-exchange. This is made explicit by our boson-pseudospin representation: ignoring the 
backscattering term in Eq. ()13|) . the Hamiltonian is identical to the spin-sector of a forward- 
scattering = KLM. Within their Abelian bosonization description, Honner and Gulacsi 
found the forward-scattering z-exchange term in the KLM directly responsible for mediating 
the double-exchange between the localized spins.— The origin of this double-exchange term 
is therefore identical to our segregating interaction. 

The shifting of the FKM's charge-sector Bose frequencies produces distortions of the c- 
electron density in response to the local /-occupation [see Eq. (fTTj) ]. These deviations from 
the homogeneous noninteracting density may be interpreted as polaronic objects;^-^ note 
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however that because of the lack of fluctuations in the / orbitals, these distortions are frozen 
into the ground state. This illustrates an important departure from the KLM physics, where 
the spin-flip {J±-) exchange terms cause the ^-component of the lattice spins to fluctuate, 
giving the distortions of the c-electron spin density (i.e. spin polarons) mobility. 

It is possible to modify the FKM in order to replicate this aspect of the KLM physics. 
The simplest such extension is an on-site hybridization term between the c and / orbitals, 
"^hyb = V J2j{^]fj + H.C.}: adding Ti^yb to Eq. (|T)) gives the quantum Falicov-Kimball model 
(QFKM). Using a bosonization mapping at the Toulouse point, Schlottmann found that the 
J_L-exchange term of the Kondo impurity is equivalent to the hybridization potential in the 
single-impurity limit of the QFKM.— In the lattice case, the polaronic distortions acquire 
mobility as in the KLM: this coupling of the c- and /-electron densities may be identified as 
a Toyozawa "electronic polaron" .— Electronic polarons in the QFKM have previously been 
studied by Liu and Ho;'^ our work on the ID QFKM largely confirms their scenario.— A 
complete account of this work is in preparation.— 

B. The longitudinal fields 

The other two terms in the transformed Hamiltonian involving the r pseudospins are a 
constant and a site-dependent longitudinal field. The former is only of importance to the 
VTP: the renormalization of the /-level by the Coulomb interaction will drive a "classical" 
valence transition. The sign of this term is proportional to np — |, which implies a strong 
dependence upon the noninteracting band structure: if > 0, the /-level will be renor- 
malized upwards, emptying its contents into the c-band; for e/ < the /-level is lowered 
below the c-band, and so all electrons will eventually possess /-electron character. Since 
this term does not determine the configuration adopted by the / electrons but rather only 
their number, we leave further discussion to when we analyze the VTP phase diagram. 

The site-dependent field is of more general interest. This originates from the 2k p- 
backscattering correction and is directly responsible for the well-known crystalline order in 
the FKM. Before demonstrating how the crystalline /-configurations can be extracted from 
this term, we first briefly review the present understanding of the weak-coupling periodic 
phases. 

The origin of the crystalline order is a Peierls-like mechanism: a one-dimensional metal 
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is always unstable towards an insulating state when in the presence of a periodic potential 
with wavevector 2kF-^^ In the context of the FKM, a Peierls instability can arise when the 
/-electrons crystallize in a periodic configuration with wavevector 2kF. This is the case 
for weak coupling and we repeat a theorem due to Freericks and Falicov:^^ given rational 
c-electron density = p/q {p prime with respect to q) and G/t ^ g, then for /-electron 
density n-^ = p-^/q {p^ not necessarily prime with respect to q) the / electrons occupy the 
sites X = nq + kj where n is an arbitrary integer and the kj satisfy the relation 

(pkj) mod q=j, j = 0, 1, . . . , - 1 (25) 

For example, consider the case ra'^ = | and n-^ = |. The unit cell has eight sites, and the / 
electrons occupy the first, second, fourth, fifth and seventh sites [see Fig. ((21) (a)]- 

Our approach reproduces this important result. In the pure crystalline phase, the c- 
electron spectrum is gapped.— We therefore replace (pixj) in the cosine term of Eq. (jlSp 
by the uniform average (0). Ignoring the Ising interaction, at weak coupling the r 

pseudospins are therefore arranged by the field 

^ J2 r! cos (2[(0) + /C(j) - k^x,]) (26) 

The string operator /C(j) is a constant of the motion and so it may be replaced by its 
eigenvalue. Referring to Eq. (fT^. this term subtracts the magnetization of the r pseudospins 
more than a to the right of site j from the magnetization of the r pseudospins more than a 
to the left of site j: for an infinite chain in the pure crystalline phase this quantity vanishes, 
/C(j) = 0. The value of (0) must be chosen to minimize the backscattering energy. This of 
course implies a non-trivial dependence upon the /-electron concentration: for a Q-site unit 
cell with /-electron concentration q-^/Q, we must minimize (0) over the sum 

l=l,qf l=gf + l,Q 

where A; = cos{2[(0) — kpXi]} and / is chosen such that Ai < A2 < . . . < Ag, lying within 
the unit cell. This minimization is most easily accomplished numerically; (0) is restricted 
to values in the interval [0, vr). 

The sum Eq. ^I7\ assumes that the q^ f electrons per unit cell will occupy the q^ lowest- 
energy sites in the potential Eq. (j26|l . In terms of the pseudospins, there is a fixed magneti- 
zation q-^ /Q — \ per unit cell; the j-spins occupy the sites with the smallest magnetic field, 
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FIG. 2: (a) The configuration of the /-ions in the weak-couphng homogeneous unit cell for n"^ 



n 



f 



. The filled and empty circles represent occupied and unoccupied sites respectively, (b) The 



pseudospin representation for the configuration (a) . (c) The variation of the effective magnetic field 
/left produced by the 2A;i?-backscattering correction across the unit cell, h^g is in units of 2GAa/a. 

with the |-spins sitting on the remaining Q — sites. For the example above with '^'^ = | 
and = |, we find (0) ~ 0.589 and so the / electrons experience a potential 



2GAa 



a 



^r/cos (^1.178 



Snxj 
4a 



We plot this potential in Fig. ©(c) along with the r pseudospin orientations [Fig. ©(b)]. 
It is in good agreement with the exact result Eq. (f^ . although there is some ambiguity 
with respect the position of the /-electron at the fifth and sixth sites. Closer correspon- 
dence may be achieved by taking into account higher-order backscattering processes; because 
bosonization is fundamentally a long-wavelength method, this approach does not replace the 
exact calculations. Nevertheless, our analysis convincingly demonstrates that bosonization 
is capable of describing crystallization of the / electrons in the FKM. 

Before proceeding to a discussion of the FKM's phase diagram, we note that as the number 
of / electrons is limited in the FKM, it may not be possible for a pure crystalline phase 
to gap the c-electron spectrum at the Fermi energy. In particular, for irrational c-electron 
filling the field Eq. (j2EI) will be incommensurate with the lattice. Although this situation 
remains unclear, in the related case of rational c-electron filling n'^ = p/q and /-electron 
filling p'/q < < {p' + l)/q the system phase separates into regions with periodic phases 
determined by Eq. for p^ = p' and p^ = p' + 1.— It is also known that for /-electron 
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concentrations < 0.371 and > 0.629 the system exists in a mixture of a crystalline 
phase and a homogeneous phase (the empty or full configuration). This phase separation 
behaviour cannot be explained purely in terms of c-electron backscattering. 



IV. THE EFFECTIVE HAMILTONIAN 

By themselves, the Ising interaction [Eq. (j^Uj) ] and the longitudinal field [Eq. (j^Uj) ] explain 
the SEG and crystalline phases respectively. To understand the origin of the phase separa- 
tion or interpret the numerically-determined phase diagram, however, we must consider the 
interplay of these terms. In particular, it is desirable to have a simple effective Hamiltonian 
for the / electrons that includes both the crystallizing and segregating tendencies of the 
FKM. 

The transformed Hamiltonian Eq. (|T8|) offers a straight-forward route to such a descrip- 
tion of the /-electrons. With the removal of the term describing the forward-scattering 
interaction, the only coupling between the two species is in the 2fcj7'-backscattering correc- 
tion. In the weak-coupling crystalline phases it is possible to completely decouple the / 
electrons from the c electrons by replacing the bosonic field by its expectation value. 

Combining Eq. (f^Bj) with the interaction Eq. ()20|) . we obtain an effective spin-i Ising model 
for the / electrons valid throughout the region where the crystalline phases are realized: 

2 9C An 

= -t: V T;5^ixj - x,,)t^, + V r/ COS (2[(0) + /C(j) - kpx,]) (28) 

This is an important result, but our approach is not limited only to the crystalline phases: 
for other configurations, the form of the effective Hamiltonian Eq. (PH|) remains valid, al- 
though the site-dependence of the longitudinal field is different. In the crystalline phases 
/C(j) is vanishing; in the segregated or phase separated states, however, /C(j) has linear 
variation. Ignoring the short-range deviation of sgn„(xj) from the true sign- function, we 
write 

m - ^ E fen - r^-n) (29) 

n=l 

Assume a phase separation between phase A and phase B with the boundary at j = 0. Then 
for j' ^ 1 we have approximately^ 

Cn 

)C{f)^m + — i{r')A-{r')B)\f\ (30) 
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where the subscripts A and B refer to the magnetization in the A and B phases respectively. 
We have chosen the sign of the hnear term by choosing phase A to be reahzed to the right 
and phase B to the left of j = 0. Note that /C(j) is constant for a pure phase, as we expect. 

For the SEG phase, the division of the lattice into empty and full sections implies a 
variation /C(j) ~ {Ga/2vp)\j\. Although the conduction electron spectrum does not display a 
gap, the decoupling procedure for the field 4>{xj) used in Sec. IIII BI mav be easily generalized. 
In the SEG phase, the conduction electrons are restricted to a fraction (1 — n-^) of the lattice, 
where they behave as a non-interacting electron gas. We therefore replace in Eq. ()26p by 
its non-interacting average (0) = to obtain the effective pseudospin Hamiltonian in the 
SEG phase. 

The phase-separation between the crystalline and empty phase is the most complex sit- 
uation to analyze, as we must account for the very different behaviour of the c-electrons for 
the two configurations. Exact diagonalization calculations on 3200 site chains reveal that 
the momentum distribution of the c-electrons is essentially a superposition of the gapped 
and noninteracting forms corresponding to the crystalline and empty regions of the lattice, 
with vanishingly small correction due to the interface of these phases in the thermodynamic 
limit — > oo.^^ Decoupling the c-electron fields as above, we take different averages of the 
0-field in the bulk of the two phases: for the empty phase, we use the noninteracting value 
{(p) = while we determine (0) for the crystalline phase as in Sec. IIIIBl Approximating 
/C(j) as in Eq. (j30|) . we have /C(j) ~ {Ga/2vF)'m\j\ where m = {t'')p + | {{t'^)p is the 
magnetization of the periodic phase). 

The effective Hamiltonian across the phase diagram is ferromagnetic Ising model in a 
oscillatory longitudinal field. This model exhibits all the most important aspects of FKM 
physics. It is, however, important to add here a cautionary note about the limitations 
of Eq. (j^Hj) . Bosonization is an inherently long- wavelength method, and so it is therefore 
unreasonable to expect Ti^ff to precisely reproduce the microscopic details of the /-electron 
configuration realized for given G, n'^ and . Rather, Ti^g is primarily relevant to the 
long-wavelength physics, with the site-dependent longitudinal field acting as an essentially 
approximate account of the short-range crystallizing interactions. Furthermore, the form of 
n is rigorously quantitatively valid only for G ^t. Nevertheless, we expect that TY^ff is at 
least qualitatively correct over a much larger region of the phase diagram.— 
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V. THE GROUND STATE PHASE DIAGRAM 
A. The crystallization problem 

In the CP, the concentration of the /-electrons is fixed: the problem of the ground state 
phase diagram is then reduced to finding the pseudospin configuration with magnetization 
n-^ — i that minimizes the energy {H^g). This lattice-gas problem, although conceptually 
simple, does not have a general solution. It is therefore appropriate to use approximate 
methods to understand the physics. 

In general, the segregating Ising interaction has a range a that extends over several lattice 
sites. To understand the segregation, however, we need consider only the nearest-neighbour 
value of the potential ^^(a). That is, we write the Ising interaction 



This is justified as for realistic cut-off A.^{k) the interaction potential Sa{x) is attractive for 
X < a but falls off very quickly with distance [see Fig. (^]. Truncating the interaction should 
not significantly alter the critical properties of the model, while considerably simplifying the 
analysis. 

The weak-coupling phase separation into the empty and a periodic configuration requires 
us to extend the Ising interaction beyond the nearest-neighbour approximation used above. 
The restriction to fixed magnetization makes this a challenging problem to analyze and a 
general criteria for the phase separation is beyond the capabilities of our approach. We 
shall nevertheless demonstrate the origin of the phase separation for a single set of input 
parameters, noting the importance of considering a delocalization length a > 2a. 

1. Segregation 

In the pseudospin "language" of the effective Ising model Eq. (j2Hl), the SEG phase cor- 
responds to a ferromagnetic state with two domains: a single block of j-spins occupying a 
fraction of the lattice and |-spins in the remaining (1 — n^)N sites. From the form of 
cff, we see that the critical Coulomb repulsion Gc for the onset of segregation is related to 




(31) 
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the critical ratio J/ h for the onset of ferromagnetism in the model 

n = -jJ2 + hY,r- cos{ujj + 0,) (32) 

j j 

where ujj and (pj take different constant values in different macroscopic regions of the lattice. 

The Ising model Eq. ()H2j) has been studied for constant u and by Sire^. For u/ir 
irrational, i.e. the quasiperiodic Ising Model (QPIM), it is found that the critical Ising 
coupling has the form = h/ sinl^uj). At couplings J > Jc the ground state is ferro- 
magnetically ordered; the adiabatic phase (where the spins align antiparallel to the di- 
rection of the longitudinal field hj) is however only realized for J < Jc2 < Jc where 
Jc2 = hsm{^r]uj) sin(|[r7 + Iju;)/ sin(|u;) and r] is the largest integer smaller than tt/lu. For 
the intermediate couplings Jc > J > Jc2 the ground state is a "mixture" of the adiabatic and 
ferromagnetic phases. This "mixed" state consists of a quasiperiodic arrangement of clus- 
ters of adiabatically- and ferromagnetically-ordered spins. These clusters form as neither the 
Ising interaction nor the magnetic field are strong enough to order the entire lattice: ferro- 
magnetic clusters occur where the magnetic field is weak compared to the Ising term, while 
paramagnetic clusters are found where the Ising term is weak compared to the magnetic 
field. Note that this is not a phase separation phenomenon. 

The work on the QPIM in Ref. |22l was performed within the grand canonical ensemble 
and so we must be cautious in relating these results to the effective Hamiltonian Eq. ()28p. 
The expression for the critical Ising coupling Jc was deduced from general arguments that 
should remain valid at fixed magnetization. Indeed, the difference in energy per site between 
the single-domain and the two-domain (SEG) solution vanishes as 0{N~^) in the thermody- 
namic limit. The QPIM results should therefore correctly capture the competition between 
the adiabatic and ferromagnetic orders present in Eq. ()28j) : this provides a condition for 
segregation to dominate crystallization. Although the QPIM at weak- and strong-coupling 
corresponds to the behaviour seen in the small- and large-G FKM, the agreement breaks 
down at intermediate coupling. This is due to the use of the grand canonical ensemble 



in Ref. 



tior^ 



371 . as phase separation cannot occur in the FKM without fixed electron concentra- 



We note that the effective Ising Hamiltonian derived for the SEG phase must always 
display ferromagnetic order. That is, within its range of applicability, the Ising interaction 
always dominates the magnetic field. Although the SEG effective Hamiltonian might dis- 
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play adiabatic order at weak coupling, since crystallization in the FKM occurs in this limit a 
different form of the longitudinal field must be used in Eq. (j2HI)- To use the QPIM condition 
to determine the boundary of the SEG phase, we therefore assume that the range of appli- 
cability of the SEG effective Hamiltonian corresponds exactly to the extent of ferromagnetic 
order. 

A further difficulty encountered when applying the FM condition derived for the QPIM 
is that the frequency of the magnetic field in the SEG phase takes two values uj± ~ 
2(7r?2^ ± Ga/2vp) for each bulk phase (i.e. the empty and full sections of the lattice). Al- 
though it is not possible to determine which value is realized for which section, the FM 
condition also holds for half-spaces and so we choose uo^ which gives the observed mono- 
tonic dependence of the critical line on the filling for weak- and intermediate-couplingi^. 
Since bosonization is quantitatively correct in the weak-coupling limit, the expressions for 
uoj^ will only rigorously hold for G small as compared to the conduction electron bandwidth. 
In this limit segregation occurs for n'^ — > and so uoj^ Gajvp] we use this form to determine 
the weak-coupling critical line. 

Comparing the coefficients in Eq. (j^H|l with those in Eq. ()H2j) we find after some algebra 
the condition for segregation 

G (XL . . , G^d J . QiAtt , . 

lim sin(c<j4./2) = sin(Gca/2t>i7') = — ^ , , (33) 

n--o Vf ^ +/ ^ \ c I t) ^^^^^^ V j 

We immediately deduce an important feature of the phase diagram: from general principles 
we know that a > Oik^^ as rf 0. Since Sa{a) ~ for a ^ a, the denominator of 
the RHS of Eq. (j^Hj) tends to a constant as a ^ cxd. For the expression to be consistent, 
we hence require Gca/vF = constant > as n'^ ^ 0: we recover the result that segregation 
occurs at arbitrarily small G in the limit of vanishing conduction electron concentration^. 

For finite filling, we can use our estimate for a/a ~ \/t/G to obtain the critical line Gc 
at weak- to intermediate-coupling. Assuming Gca/vp ^ 7r/2 for small G^ we linearize the 
sine function in Eq. (jH^ : after some algebra we find 



G, = 4tsm{7Tn')^l^^ (34) 

At weak coupling, we have a ^ a: Taylor-expanding the RHS of Eq. (jH^ in powers of a/a, 
we keep terms up to second order. The coefficients in this expansion are dependent upon 
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FIG. 3: (color online) Dependence of Gcd/vp on Gc/t for three values of the ratio r. The data is 



taken from Ref. 
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the form of A„(/c) used; for exponential cut-off we have 



2a 



(35) 



Substituting our estimate for a into this equation, we thus expect a hnear relationship 
between Gcd/vp and Gc- This also holds at intermediate couplings, as clearly verified by 
the numerical results of Ref. Q: we plot Gca/vp as a function of Gc in Fig. Q for three 
values of the fraction of electrons in the c-band r = n'^/n = 0.75, 0.5 and 0.25. After some 
re-arrangement of Eq. (j35|) . we obtain the general form of the critical line 

G.(r.n)^ /^':;7'7"' (36) 

1 — 2C [rj sm[TTrn) 

The numerical constants B{r) and G{r) are the intercept and gradient of the lines 
in Fig. (jni) respectively; they are related to the fitting parameters in Eq. hj A = 
B{rf/47T^ and a/a = ^ B{r)t/2C{r)G^. 
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FIG. 4: (color online) Dependence of Gc/t on n for three values of the ratio r. The data is taken 
from Ref. Iisl. The critical lines Gc{n) of best fit are as derived from Fig. (jS)). For each value of r, 
the SEG phase occurs for G > Gc{n). 

From our fit to the lines in Fig. Q we find the critical lines for the three values of r 

0.5666 sin(n7r/4) 



1 - 1.5212 sin(n7r/4) 

0.5sin(n7r/2) 
1 — sin(ra7r/2) 

0.4112 sin(3r27r/4) 



0.25 
0.5 
0.75 



(37) 



U- 0.8982 sin(3n7r/4) 

These three curves, along with the numerical data, are plotted in Fig. (0)). The curves 
track the data very well for both r = 0.25 and r = 0.5; for r = 0.75, however, there is a 
significant divergence between Eq. (j37p and the numerical results as the coupling increases. 
The curve Eq. fj37|) has a maximum at n = 2/3 (i.e. = 1/2), but no evidence of this 
maximum is found in the weak-coupling numerical results. Rather, we expect the critical 
line to continue to diverge as half-filling is approached. We thus conclude that there is a 
change in the form of Gc{n) at intermediate-coupling. The numerical analysis of Gruber 
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et al. indicates that this occurs at approximately G ~ 2.5^^^, which is consistent with the 
observed deviation from the weak-couphng critical line in Fig. Q. 

Any deviation from the weak-coupling form Eq. (j36|) for r = 0.25 and r = 0.5 is much less 
obvious. Since the G ^ oo asymptotic form of Gc{r, n) stated in Ref. Q is not the same as 
that given by Eq. (j37|) . we do however expect that a different expression is valid aX G ^ t. 
A new functional dependence on n in the strong-coupling regime is reasonable and does not 
contradict our own analysis: we have emphasized that bosonization is only quantitatively 
accurate for weak-coupling. Importantly, the physical processes driving the segregation will 
remain invariant across the phase diagram. 



2. Phase separation 

At weak-coupling and sufficiently small or large /-electron concentration, the FKM is un- 
stable towards a phase separation between a homogeneous and a crystalline configuration.— 
Unlike the SEG phase, it is necessary to consider the range of the forward-scattering Ising 
interaction as extending beyond the nearest neighbour to observe these phases. 

To illustrate the importance of these long-range terms, we examine the weak-coupling 
limit of Eq. ()28|) for the case n'^ = |, n-'^ = i. Of the Ising interaction Eq. (j20|) we keep 
the nearest-neighbour Ji and next-nearest neighbour J2 terms. For G ^ t, we discard 
/C(j) in the cosine's argument, leaving a staggered-field variation. We thus find an effective 
Hamiltonian of the form 

^ = -^1 E - ^2 E - h J2i-iyr^ (38) 

j j 3 

We calculate E = {Ti) for three situations: (a) the most homogeneous phase with period-4 
pseudospin unit cell [Tiii]; (b) phase separation between the empty phase {[[]) and the 
period-2 phase with unit cell [][]] and (c) segregation. We find the energy per site for each 
of these configurations 

f 

— \h config. (a) 

E/N = \ -1 - Ih config. (b) (39) 

^-\Ji-\j2 config. (c) 
These expressions hold in the thermodynamic limit N ^ 00. 
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We see that even an arbitrarily small J2 destabilizes configuration (a) toward phase 
separation. Since at weak-coupling we expect J\ <^ segregation will however not occur; 
instead configuration (b) is the most stable. This is the weak-coupling phase separation 
between a crystalline and the empty phase found by Freericks et alP Although our analysis 
does not provide a general condition for this peculiar form of phase separation to occur, 
it shows that the physical origin of this effect is the competition between segregation and 
crystallization. We also see how fixed c- and /-electron populations are essential for the 
appearance of this phenomenon. 

The phase separation between periodic and uniform states in the n < 1 FKM is not con- 
fined to the weak-coupling limit of the phase diagram, but is also present at intermediate- 
and strong-coupling . ^^1^^ For G > t our bosonization approach will at least qualitatively 
capture the physics of the FKM: we therefore expect that the competition between segre- 
gation and crystallization that we have identified as the origin of the weak-coupling phase 
separation will also be responsible for these intermediate- and strong-coupling phases. 

B. The valence transition problem 

In contrast to the CP, the VTP has received little attention in the FKM literature, 
despite the two interpretations being very closely related. In both the CP and the VTP 
the /-orbital occupation is a good quantum number, and the ground state may be defined 
as the configuration of the / electrons that minimizes the energy of the c electrons. The 
only difference between the CP and the VTP is that in the former the distribution of 
the electrons across the orbitals is fixed, while in the latter the interactions determine the 
equilibrium populations. 

For given interacting equilibrium populations in the VTP, the configuration adopted by 
the /-electrons should be the same as in the CP for the same fixed electron populations. As 
discussed in Sec. IIII Bl the first two terms of Eq. ()18|) determine the equilibrium distribution 
of the electrons across the c and / orbitals. They can be identified as the noninteracting 
c-electron Hamiltonian and the /-level shift due to the Coulomb repulsion. To estimate 
the electron distribution for finite G we therefore assume that the distribution of the nN 
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electrons across the two orbitals in the FKM is the same as in the system 



<;,+i + H.c.| + [e, + G(n5 



(40) 



3 



j 



That is, the contribution of the ordering terms in Eq. ()18|) to the shift in electron density 
between the c and / orbitals is taken to be negligible. This can be easily justified for a 
thermodynamically large system: the difference between the energy per site of the SEG 
phase and the empty or full phases due to the Ising interaction is of order and the 

average value of the bacscattering longitudinal field across the lattice is vanishing. 

We find that for noninteracting c-electron population (fixed by the band structure) 
the c-electron population in the interacting system is given by 



In Ref. phase diagrams for the CP in the n^-n-' plane at constant G are presented For 
each G the boundary between the SEG phase and the crystalline or phase separated states 
is given by a straight line of the form n'^ = 7(1 — n-^) where 7 is a constant determined 
from the numerical phase diagrams. Using the fixed electron concentration condition we 
may re-write this 

n' = ^— (1-n) (42) 
1-7 

where n = + n-^ is the total electron concentration. Equating the RHS of Eq. fl41|) 
and Eq. we obtain the equation 



where tq = rig/n is the fraction of c-electrons in the non-interacting system. For given ro 
and G, the value of n that solves Eq. ()43p gives the maximum filling for which the SEG 
phase is stable. 

using t..s p.oced.e we calculate .o,„ the results of Ref. U the cnt.cal value of „ fo. 

tq = 0.25, 0.5 and 0.75. Proceeding with our analysis as in the CP (Sec. IV A 1|) we find that 
the linear relationship between Gcd/vp and Gc is also well obeyed in the VTP [Fig. (jSJ]. 
As such, at weak-coupling the critical line Gciro,"!^) has the form given by Eq. ()36p. In 




(41) 




(43) 
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FIG. 5: (color online) Dependence of Gca/vp on Gc/t for three values of the ratio rg in the VTP. 
The data is taken from Ref. 2. 

particular, we find from the linear best fit to the data in Fig. (jSj) the following expressions 

0.389 sin(ra7r/4) 



Gc{n)/t 



1 - 1.289 sin(n7r/4) 

0.363 sin(n7r/2) 
1 - 0.992 sin(n7r/2) 

0.307 sin(3ri7r/4) 



ro = 0.25 



ro = 0.5 



ro = 0.75 



(44) 



^ 1 - 0.923 sin(3n7r/4) 

These are illustrated in Fig. (jHl) along with the critical lines for the associated r = tq CP 
[Eq. (jHTj) ]. As before, we find excellent agreement between the data points and the fitted 
curves for both tq = 0.5 and tq = 0.25. Again, however, we find for tq = 0.25 a significant 
divergence between the curve Eq. and the data for higher values of G. The origin of 
this discrepancy is presumably the same as in the CP. Note also that only four numerical 
values are presented for ro = 0.25: for G > ?>t the segregated configuration is realized for 
all < 1. This is not unexpected, as in this case we have the smallest np, and thus largest 
shifting of the /-level, for given n. 

As illustrated in Fig. the division of electrons between the two orbitals in the CP 
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FIG. 6: (color online) Dependence of Gdt on n for three values of the ratio r for the VTP. The 
data is taken from Ref. I la. The critical lines Gc{n) of best fit (thick dashed lines) are as derived 



from Fig. ©. For each value of tq, the SEG phase occurs for G > Gc{n). The thin dotted lines 
are the critical lines in the r = tq CP. 

strongly affects the position of the critical line for segregation: the more / electrons relative 
to c electrons, the smaller the value of G required to cause segregation. For < |, turning 
on the interaction in the VTP will shift the /-level to a lower energy relative to the c-electron 
band, thus causing a transfer of electrons from the c to the / orbitals. Accordingly, we find 
that segregation in the VTP occurs at a lower value of G than in the r = rg CP [Fig. ©]. 
Eventually, the /-level will be shifted below the bottom of the c-electron band; this is the 
case for couplings 

G,.„ > (45) 

The absence of any c electrons to cause crystallization or segregation means that any /- 
electron configuration is the ground state. For tq = 0.25 and n = 1 the critical coupling 
Gfuii ~ 2.34t, explaining the absence of any G > 2t data. 

Conversely, for > and hence > |, turning on the interaction will shift the /-level 
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to higher energies and empty the /-electron orbitals. Segregation may not occur in this case 
at all, and the large-G configuration is the empty phase. From our analysis, we estimate 
that this will be realized for 

-4t[cos(n7r) - cos(ron7r)] . 
Ge.p.. > (46) 

Note that we assume n> \. This scenario is strongly supported by Farkasovsky's numerical 
study of the n = 1 VTP.— In his ef-G phase diagram, he found that for e/ > (ng > |) all 
electrons occupy the c-orbital states for sufficiently large G, while for e/ < < i) the / 
orbitals become fully occupied as G is increased. 

We note in concluding that we have not addressed the case where the /-level does not 
lie at the Fermi energy in the noninteracting system. For example, for e/ = and n < ]^ 
the /-level will lie a finite energy above ep = —2tcos{n7c). On the basis of our analysis, it 
appears that a non-zero /-population will eventually appear as the /-level is shifted in the 
presence of a finite G. As G is further increased, the c-electron band is eventually emptied 
into the /-level. Turning on the interactions, we thus evolve from a state without any / 
electrons into a state with all electrons in the / orbitals. We must regard this result with 
caution: since bosonization is an effective field theory for the excitations about the Fermi 
energy, it is difficult to include the localized electrons whenever ej ^ ep- In particular, the 
bosonic wavelength limit a implies an effective bandwidth cut-off for the excitations about 
cp- What happens when e/ lies outside of this effective bandwidth is unclear and we must 
go beyond the framework of bosonization to understand this situation. For this reason, from 
the point of view of bosonization the VTP is a more challenging problem than the CP. 



VI. EXTENSIONS OF THE FKM 



The FKM is often studied in a modified form with the addition of extra terms to the 
basic Hamiltonian Eq. (0). The most common extensions are c-f hybridization , ^i^i^i^^i^^ /- 
electron hopping^S*^ or the introduction of spin.— i^Siii In the first two cases, the extension 
has a dramatic effect upon the physics: the occupation of each localized orbital is no longer 
a good quantum number. The CP and VTP results are applicable only as limiting behaviour 
and we cannot easily incorporate these additional terms into the analysis presented above. 
As such, below we will briefly consider several extensions that maintain the "classical" nature 
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of the / electrons: intraorbital nearest-neighbour interactions and the addition of spin. Our 
bosonization formalism is very well suited to assessing the impact of these extensions upon 
the ground states of the "bare" FKM. 



A. Nearest-neighbour interactions 

Our study of nearest-neighbour interactions is confined to their effect upon the CP results. 
The same conclusions also hold for the VTP so long as the densities are normal-ordered. 

1. c electrons 

We write the nearest-neighbour interaction between the c electrons 

Kc = K5;]n^SVi (47) 
j 

It is sufficient here to examine only the forward-scattering contributions of this interaction 
as the Umklapp and backscattering contributions are only relevant at half-filling^^. Since 
we are only interested in the effect of Ticc on the long-wavelength physics of the FKM, we 
may apply the continuum-limit approximation and absorb the interaction into a free-Boson 
Hamiltonian 

^0 = E { (^^^(^.))' + {dj{x,)y\ (48) 

j 

The new Bose fields are related to the K = fields by the relations 



0(x,) = -j=<p{x,), e{xj) = VKe{xj) (49) 
V K 

(50) 



K 



where 



K = ^ = (51) 

y/1 + 2VcVF/{nat^) 

The details of this rescaling procedure are identical to the argument for the forward- 
scattering sector of the XXZ chain^. Note that for an attractive interaction Vc = —Ttat^/vp 
the velocity of the Boson modes vanishes: this indicates the break-down of the bosonization 
method, as the Luttinger liquid is unstable towards the clustering of the c-electrons. We 
may expect the SEG phase to be realized whenever this condition holds. 
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The bosonized FKM with the interaction Eq. ()47|) is identical in appearance to the = 
bosonized Hamiltonian [Eq. (fT^ ]: the first term in Eq. ()13|) is however replaced by Eq. (jlH)) 
and the Bose fields in the other terms are replaced by their scaled forms Eq. (|49|) . Our 
analysis of T^fkm + T^cc also proceeds in a similar way to that in Sec. irm although we rotate 
the Hilbert space using a different canonical transform 

f/ = exp|.^X;r;^(x,o| (52) 

As before, we find the effective segregating interaction Eq. (j^Uj) . but with the coefficient 
changed by a multiplicative factor 

G^a" Co" 1 , , 

(53) 



2'kvf 2ttvp 1 + 2Vca/TTVp 

As we can see, the effect of a repulsive (attractive) nearest-neighbour interaction between 
the c electrons is to suppress (enhance) the segregating interaction. This conclusion is not 
surprising: the interaction Eq. (jTfj) rescales the charge compressibility n of the c electrons 

K = Kq = KqK 

V 

where kq is the compressibility for Vc = 0. Repulsive interactions {K < 1) reduce the 
compressibility and vice versa. In the SEG phase the density of the c electrons is enhanced 
due to their confinement to a fraction 1 — n-^ of the lattice. As such, a reduced (enhanced) 
c-electron compressibility will resist (assist) the formation of this state. 

We can easily judge the effect of Hcc on the position of the critical line Gc- Following the 
same arguments as in Sec. IV A If we find in the limit n'^ the asymptotic form 



a,^(^±^'-}^).J^^ (54) 

The term under the square-root is constant for small c-electron fillings; the expression in 
brackets therefore determines the small-n'^ form of the critical line. We thus find that for 
Vc > 0, the SEG phase is only realized above a finite Coulomb repulsion even in the limit 
of vanishing c-electron concentration. For attractive interactions, however, the system is 
unstable towards segregation for any c-electron filling such that K^^ = 0. 
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2. f electrons 



The nearest-neighbour interaction between the / electrons is much easier to analyze. It 
is a simple matter to write the interaction term in the pseudospin representation 

^// = VfJ2 = ^/ E ^N+i (55) 

j j 

This nearest-neighbour Ising interaction may be immediately incorporated into our effective 
pseudospin model Eq. (j2Hl)- 

Quite clearly, an attractive interaction potential Vf < will make the system unstable 
towards the SEG phase even for G = 0. Crystallization may still occur, although only at 
finite coupling strength. Furthermore, we expect that the SEG phase will be realized even 
at half-filling for sufficiently large G. We consider a large-G expansion where we project the 
FKM into a truncated basis excluding simultaneous occupation of both the c and / orbitals. 
To first order in G^^ we find the effective strong-coupling Hamiltonian 



{^+yf)j:-Fu^ (56) 



where = |(n^ — nj) and the magnetization is fixed at = |(n^ — n^). The sign of the 
nearest-neighbour interaction is ferromagnetic for G/t > —2t/Vf implying the formation of 
the SEG phase. Of course, higher-order [at least (9(G'^^)] terms complicate this analysis, 
but by increasing G we can make their contribution arbitrarily small. 

More interesting is the case of a repulsive potential V/ > 0. Here Eq. ()55|1 hinders 
segregation, and for sufficiently strong Vf may suppress it entirely. This is dependent upon 
the sign of the nearest-neighbour Ising interaction in the TYfkm + ^// effective pseudospin 
model: the SEG phase cannot be realized unless the nearest-neighbour Ising interaction is 
ferromagnetic. By equating Eq. (jHT|) and Eq. ()55|) we immediately find a condition for the 
appearance of segregation: 

Since for n'^ <^ 1 we have Sa{ci) ~ ^ 0{kp), the RHS of Eq. ^7\i should be finite at low 
c-electron filling. Phase separation as discussed in Sec. IV A"2] will nevertheless still occur as 
the range of the interaction Eq. ()20j) extends beyond nearest-neighbours, with these higher- 
order terms in the pseudospin model remaining unaffected by the addition of The 
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most important of these extra terms is the next-nearest-neighbour interaction: if the condi- 
tion Eq. (j37j) is not satisfied, this term is the dominant ferromagnetic couphng, and hence 
orders the /-electrons into a single cluster where only every second site is occupied. That is, 
we may expect that a large portion of the SEG phase in the phase diagram Fig. (@)) will be 
replaced by a phase-separation between the empty and period-2 crystalline configurations. 

Numerical results for the FKM with nearest-neighbour /-electron repulsion confirm this 
scenario: Gajek and Lemahski have studied the effect of Eq. ()55|) in the canonical ensemble 
for Vf = O-IC^^ For a repulsive potential of this form, the SEG phase was not realized at 
any coupling strength or electron filling. For -C 1, the /-electrons indeed phase separate 
into the period-2 crystalline and empty configurations. Interestingly, for the Uq = case 
presented, phase separation is realized only for n < 0.4. This indicates is a significant 
truncation of the range of the segregating interaction with increasing filling. 

B. Spin 

To use the FKM as a model of any realistic condensed-matter system, we are required 
to relax the assumption of spinless electrons. Simply adding a spin index to the fermionic 
operators in Eq. is, however, not enough: we must take into account the orbital struc- 
ture of the localized states. Because of the small radius of the / orbitals, the intra-ionic 
correlations are very strong, prompting us to introduce a Coulomb repulsion U between the 
/ electrons in our spinful model. We thus write 

^ = E E {s>S-t-i,. + H.C. } + uY: r^,,Ai + ^ E E 4.4.' (58) 

3 . 3 3 o-,.' 

We consider here only the limit U = oo where double occupation of an /-orbital is excluded 
from the physical subspace. This situation has been numerically studied by Farkasovsky 
for the fixed total electron concentration n = 1 in both the CP and VTP interpretations.— 
Since double occupation is forbidden, we may represent the / operators in terms of spinless 
fermion (ho Ion) operators 

E<. = (1-SS) (59) 
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That is, at any site without an /-electron we find a spinless hole. We hence re-write Eq. fl58|) 



+ H.c.} + G 5: 5^(1 - e]e,.)n,V (60) 



The /-orbital occupation is thus described by spinless fermions as in the usual FKM; the 
condition for fixed total electron concentration is however written n = (l/N) Ylji^ " + 
^^(n^^)]. The spin-modes of the c electrons cannot be removed as for the / electrons. 

The bosonization procedure outlined in Sec. |n] requires little modification to include the 
spin degrees of freedom. We define boson fields in terms of the density fiuctuations Pu,a{k) 
in each spin-channel 



- ^^^) = -'1^1^ ^P'^Amaik)e-''^^ (61) 

u k=iO 

= ^ E E ^^P'^Ak)K{k)e-''^^ (62) 

Boson fields with different spin-indices commute; fields with the same spin-indices obey the 
commutation relations Eq. ((Tj) and Eq. (jH)). It is convenient to split the bosonic representa- 
tion into charge- and spin-sectors, defined respectively by the linear combinations 

<Pc{,Xj) = (63) 

M^j) = -^Ihi^j) - (Pli^j)] (64) 

and similarly for the ^-fields. This will considerable simplify the bosonic representation of 
our Hamiltonian. After some algebra we find the bosonic representation for the electron 
density operator 

= tiq - ^^^^d^(j)c{,Xj) + cos[v^(/)s(x,)] cos[v^0c(a;,) - 2kFXj\ (65) 
^ — ^ TT a 

where kp = im'^ /2a. Note that the forward-scattering contribution (second term on RHS) 
is very similar as in the spinless case Eq. (jllj) : the backscattering contribution (third term 
on RHS) however involves both the spin- and charge-sector fields. 

Again we adopt a pseudospin representation for the /-orbital occupation: we define 
Tj=\ — e]ej, and so as before spin-| corresponds to an occupied orbital and spin-| to an 
empty site. Following the same basic manipulations as for the spinless case, we obtain the 
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bosonized Hamiltonian 



27r 

TT 



i 

cos[v^0s(xj)] cos["\/20c(a^j) — 2A;i?Xj] (66) 



Excluding the last term, Eq. ()66|) is identical to its spinless equivalent. Importantly, the 
forward-scattering interaction (second last term) remains in the same form as before: by 
simply shifting the charge-sector boson frequencies we may remove this term. This requires 
us to apply the canonical transform 

[/ = exp <( i^^^ > ^ TW,{xy) \ (67) 




0' 

After some algebra, we obtain the transformed Hamiltonian 

U^niJ = ^ E E {(^^'^«(^.))' + i.dM^3)f] + G{nl - i) 5^r/ 
Co" 



TTVp 



AGAa ^ cos[V2(j)s{xj)] cos[V2(j)c{xj) - 2/C(j) - 2kFXj] (68) 



a 



where /C(j) is defined as in Eq. (fTU|) . Since the c-electron spin modes do not contribute to 
the physics, we may replace 0^ by its noninteracting expectation value, i.e 0^ = (0s) = 0. 
Substituting this into Eq. (j66|l we obtain the same effective Hamiltonian as for the spinless 
FKM. This allows us to draw an important conclusion: for the spinful model Eq. (jHn|) with 
/- and c-electron concentrations n-^ and n'^ respectively, the configuration adopted by the 
/ electrons is identical to that adopted by the spinless / electrons in Eq. with /- and 
c-electron concentrations and respectively. This explains the appearance of phase 
separation and segregation in the numerical study of Eq. (jUnj) at n-^ + = 1.— 

The obvious extension of Eq. (|58p would be the inclusion of a Kondo-like exchange be- 
tween the c and / electrons on each site j. This model may be useful for understanding the 
properties of the manganites, which are known to display a phase separated state.-^ Such 
a model could display an interesting coexistence of spin- and charge-order; this problem 
remains to be fully addressed.—"^ 
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VII. CONCLUSIONS 



In this paper we have presented a novel approach to the study of charge order in the 
FKM below half-filling. We used a bosonization method that accounted for the non-bosonic 
density fiuctuations of the c electrons below a certain length-scale a > a; we identified a as 
characterizing the delocalization of the c electrons. This delocalization of the c electrons over 
several lattice sites favours empty underlying / orbitals in order to minimize the interorbital 
Coulomb repulsion. We demonstrated in Sec. IIII Al how this directly leads to effective at- 
tractive interactions between the / electrons, and hence the SEG phase. Using a canonical 
transform, we obtained an explicit form for the segregating interaction Eq. (|2(Jj) . Since the 
canonical transform was carried out to infinite order, this interaction is non-perturbative. 
The canonical transform is a generalization of the transform used in Schotte and Schotte's 
solution of the XEP.— We argued in Sec. IIII A II for a parallel between Eq. ()20j) and the 
double-exchange interaction in the KLM, based upon the similar orthogonality catastrophe 
physics in the single-impurity limit of both models. 

The canonical transform permitted a decoupling of the c and / electrons, yielding an effec- 
tive Ising model [Eq. (j2HI)] for the configuration of the / electrons, based upon a pseudospin- 
\ representation for the occupation of the localized orbitals. This effective model, TYch, 
clearly revealed the competition between the backscattering crystallization and the forward- 
scattering segregation. T-^^g correctly predicted the structure of the CP phase diagram: we 
obtained an expression for the critical coupling required for segregation which is in good 
agreement with the numerical results. We also demonstrated that the effective model could 
successfully account for the instability towards phase separation between a crystalline and 
the empty phase in the weak-coupling FKM. Our approach was not limited to the CP, and 
in Sec. IV Bl we considered the phase diagram of the VTP. We found that the Coulomb re- 
pulsion shifted the bare /-level, causing a "classical" valence transition. The sign of the 
/-level shift is highly dependent upon the band structure. Finally, we discussed the im- 
pact of intraorbital nearest-neighbour interactions (Sec. IVI A|l and the introduction of spin 
(Sec. IVIB|) on the charge order found in the spinless model. 

Prospects for future work are promising. We have already outlined an application of our 
method to the nontrivial extension of Eq. (Q) by the addition of an on-site hybridization 
potential, the so-called Quantum FKM (QFKM).— The crystallization is heavily suppressed 
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in the QFKM, as the dominant feature of the c-electron behaviour at weak-couphng is the 
resonant scattering off the / orbitals (mixed- valence). In contrast, segregation occurs in 
the QFKM, as the responsible orthogonality catastrophe physics remains intact with the 
introduction of the hybridization fSec. 1111 A Here, however, we expect dynamic charge- 
screening processes (in analogy to the spin-screening in the KLM), with important and 
interesting consequences that we will fully explore in a forthcoming publication.— 
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